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A partial wave analysis of the pp mass-threshold enhancement in the reaction J/ip — > ^ypp is 
used to determine: its J PC quantum numbers to be _ + ; its peak mass to be below threshold at 
M = 1832±5 9 (stat.)il? (syst.) ± 19 (model) MeV/c 2 ; and its total width to be Y < 76 MeV/c 2 
at the 90% C.L. The product branching ratio is measured to be B(J/ip — > r yX(pp))B(X(pp) — ¥ 
pp) = (9.0ti'i (stat.)^ s ' (syst.) ±2.3 (model)) x 10 -5 . A similar analysis performed on ip' — s- 'ypp 
decays shows, for the first time, the presence of a corresponding enhancement with a production 
rate relative to that for J/ip decays of R = (5.08to.l5 (stat.)^;® 7 (syst.) ± 0.12 (model))%. 

PACS numbers: 12.39.Mk, 12.40. Yx, 13.20.Gd, 13.75.Cs 



An anomalously strong pp mass threshold enhance- 
ment was first observed by the BESII experiment in 
the radiative decay process J/ip — > jpp [l[ and was re- 
cently confirmed by the BESIII and CLEO-c experi- 
ments. Curiously, no apparent corresponding structures 
were seen in near-threshold pp cross section measure- 
ments, in JS-meson decays [3|, in radiative ip 1 or T — > jpp 
decays Q, or in J/ip — > ujpp decays @. These non- 
observations disfavor the mass-threshold enhancement 
attribution to the effects of pp final state interactions 
(FSI) HH. 

A number of theoretical speculations have been pro- 
posed to interpret the nature of this structure (61-TlOj] - 
Among them, one intriguing suggestion is that it is due 
to a pp bound state, sometimes called baryonium [Icj . 
an object with a long history and the subject of many 
experimental searches [Til ]. The observation of the pp 
mass threshold enhancement also stimulated an experi- 
mental analysis of J/ip —¥ 77r + 7r~?/ decays, in which a 
7r + n~r]' resonance, the X(1835), was first observed by the 
BESII experiment [l2| and recently confirmed with high 
statistical significance by the BESIII experiment [13| . 



Whether or not the pp mass threshold enhancement and 
the X(1835) are related to the same source still needs fur- 
ther study; among these, spin-parity determinations and 
precise measurements of the masses, widths and branch- 
ing ratios are especially important. 

In this letter, we report the first partial wave analy- 
sis (PWA) of the pp mass threshold structure produced 
via the decays of J/ip — > jpp and ip' — > jpp. Data 
samples containing (225.2 ± 2.8) x 10 6 J/tp events and 
(106 ± 4) x 10 6 Tp' events fill a ccumulated in the Bei- 
jing Spectrometer (BESIII) 15] located at the Beijing 
Electron-Positron Collider (BEPCII) (HI are used. 

The cylindrical core of the BESIII detector consists 
of a helium-gas-based drift chamber (MDC), a plastic 
scintillator Time-of-Flight system (TOF), and a CsI(Tl) 
Electromagnetic Calorimeter (EMC), all enclosed in a 
superconducting solenoidal magnet that provides a 1.0-T 
magnetic field. The solenoid is supported by an octagonal 
flux-return yoke with resistive plate counter muon iden- 
tifier modules (MU) interleaved with steel plates. The 
solid angle for the charged particle and photon accep- 
tance is 93% of 47T, and the charged particle momentum 
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FIG. 1: The pp invariant mass spectrum for the selected 
J/xp "/pp candidate events, (a) The pp invariant mass spec- 
trum; the open histogram is data and the dashed line is from 
J/ip — > "jpp phase-space MC events(with arbitrary normaliza- 
tion), (b) An M 2 (yp) (horizontal) versus M 2 (-yp) (vertical) 
Dalitz plot for the selected events. 



and photon energy resolutions at 1 GeV are 0.5% and 
2.5%, respectively. The time resolution of TOF is 80 ps 
in the barrel and 110 ps in the endcaps, and the dE/dx 
resolution is 6%. 

Charged-particle tracks in the polar angle range 
| cos 6*| < 0.93 are reconstructed from hits in the MDC. 
The TOF and dE/dx information are combined to form 
particle identification confidence levels for the 7r, K and p 
hypotheses; the particle type with the highest confidence 
level is assigned to each track. Photon candidates are 
required to have an energy deposit of at least 25 MeV in 
the barrel EMC (| cos6>| < 0.8) and 50 MeV in the end- 
cap EMCs (0.86 < |cos0| < 0.92), and be isolated from 
antiprotons by more than 30°. 

Candidate J '/tp — > jpp events are required to have at 
least one photon and two charged tracks identified as a 
proton and an antiproton. Requirements of |£/ m i SS | < 
0.05 GeV, where U miss = (E miss - \P m i ss \), and P t 2 7 < 
0.0005 (GeV/c) 2 , where P 2 7 = 4|P miss | 2 sin 2 7 /2, are 
imposed to suppress backgrounds from multi-photon 
events. Here E m i SS and P m iss are, respectively, the miss- 
ing energy and momentum of all charged particles, and 7 
is the angle between the missing momentum and the pho- 
ton direction. A four-constraint (4C) energy-momentum 
conservation kinematic fit is performed to the jpp hy- 
pothesis. For events with more than one photon can- 
didates, the combination with the minimum x 2 is used; 
X 2 < 20 is also required. Since there are differences in 
detection efficiency between data and Monte Carlo (MC) 
simulated low-momentum tracks, we reject events con- 
taining any tracks with momentum below 0.3 GeV/c. 

The pp mass spectrum for events that satisfy all of the 
above- listed criteria is shown in Fig. HJa). There is a 
clear signal of r/ c , a broad enhancement around M p p ~ 
2.1 GeV/c 2 , and a prominent and narrow low-mass peak 
at the pp mass threshold.consistent with that reported 
by BESII [l| and BESIII [2]. The Dalitz plot for selected 
events is shown in Fig. [Tfb). 

Potential background processes are studied with an 
inclusive MC sample of 2 x 10 8 J /if) events generated 



according to the Lund model [l7|. None of the back- 
ground sources produces an enhancement at the pp mass 
threshold region. The dominant background is from 
J/ip — > ir°pp events, with asymmetric 7r° — > 77 decays 
where one of the photons has most of the ir° energy. An 
exclusive MC sample, generated according to the PWA 
results of J/ip — > n°pp at BESII indicates that the 
level of this background in the selected data sample with 
M pP < 2.2 GeV/c 2 is 3.7% of the total. The J/ip -> n°pp 
decay channel is also studied with data, and there is no 
evidence of a pp mass threshold enhancement, which pro- 
vides further evidence that the enhancement observed in 
J/?/ 1 decays is not from background. 

A PWA of the events with M pP < 2.2 GeV/c 2 is per- 
formed to focus on determining the parameters of the 
pp mass threshold structure, which we denote as X(pp). 
The maximum likelihood method applied in the fit uses 
a likelihood function that is constructed from "/pp signal 
amplitudes described by the relativistic covariant tensor 
amplitude method and MC efficiencies. The back- 
ground contribution from the ir°pp process is removed 
by subtracting the log-likelihood values of background 
events from that of data, since the log-likelihood value of 
data is the sum of the log-likelihood values of signal and 
background events [20]. Here, the background events are 
estimated by the MC sample of J/ip — > ir°pp decays de- 
scribed above. We include the effect of FSI in the PWA 
fit using the Julich formulation [||. 

Four components, the X(pp), / 2 (1910), / (2100) and 
++ phase space (PS) are included in the PWA fit. The 
intermediate resonances are described by Breit-Wigner 
(BW) propagators, and the parameters of the /2(1910) 
and /o(2100) are fixed at PDG values. In the optimal 
PWA fit, the X(pp) is assigned to be a 0~ + state. The 
statistical significance of the X(pp) component of the 
fit is much larger than 30<r; those for the other com- 
ponents are larger than 5er, where the statistical signifi- 
cance is determined from the changes of likelihood value 
and degrees of freedom in the PWA fits with and with- 
out the signal hypotheses. The mass, width and prod- 
uct branching ratio (BR) of the X(pp) are measured to 
be: M = 1832l^ 9 MeV/c 2 , T = 13 ± 39 MeV/c 2 and 
B(J/ip -/X)B{X pp) = (9.0±S;t) x 10~ 5 , respec- 
tively, where the errors are statistical only. Figure [2] 
shows comparisons of the mass and angular distributions 
between the data and the PWA fit projections. For the 
spin-parity determination of the X(pp), the ^ assign- 
ment fit is better than that for ++ or other J PC assign- 
ments with statistical significances that are larger than 
6.8ct. 

Variations of the fit included replacing the /o(2100) 
with the / 2 (2150), the / 2 (1910) with the / 2 (1950), and 
replacing both components simultaneously; changing the 
J PC of the PS contribution, as well as consideration of 
the parameter uncertainties of the /o(2100) and /2(1910), 
were performed, and it is found the changes of the log- 
likelihood values and the parameters of the X(pp) are 
quite small, except that when replacing ++ PS with ^ 
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FIG. 2: Comparisons between data and PWA fit projection: 
(a) the pp invariant mass; (b)-(d) the polar angle # 7 of the 
radiative photon in the J/ip center of mass system, the polar 
angle 9 V and the azimuthal angle <p v of the proton in the 
pp center of mass system with M p p — 2m p < 50 MeV/c 2 , 
respectively. Here, the black dots with error bars are data, 
the solid histograms show the PWA total projection, and the 
dashed , dotted , dash-dotted and dash-dot-dotted lines show 
the contributions of the X(pp), ++ phase space, /o(2100) 
and /2(1910), respectively. 



PS the event number of the X(pp) decreases by 52%. We 
also tried fits that include other possible resonances listed 
in the PDG table HJ [7? 2 (1870), / 2 (2010), / 2 (1950), 
/ 2 (2150), /j(2220), t?(2225), M23Q0), / 2 (2340) etc.] as 
well as X(2120) and AT(2370) 0, and different J PC PS 
contributions. The statistical significances of these addi- 
tional resonances are lower than 3cr. All of the parameter 
changes that are found in these alternative fits are con- 
sidered as sources of systematic uncertainties. 

For systematic errors on the mass and width of the 
X(pp), in addition to those discussed above, we in- 
clude uncertainties from different fit ranges of M pp < 
2.15 GcV/c 2 and M pP < 2.25 GeV/c 2 , different parame- 
terizations for the BW formula, as well as different back- 
ground levels. For the systematic errors of the BR mea- 
surement, there are additional uncertainties from the ef- 
ficiencies of charged track detection, photon detection 
and particle identification, kinematic fit and the total 
number of J /if) events. The total systematic errors on 
the mass and width of the X(pp) are MeV/c 2 and 
i 13 MeV/c 2 , respectively, and the corresponding relative 



systematic error on the product BR is 
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Various FSI models |6j-|8( have been proposed to inter- 
pret the pp mass threshold enhancement. Among them, a 
BW function times a one-pion-exchange FSI factor 8] can 



FIG. 3: (a) The pp invariant mass spectrum for the se- 
lected ip' — ¥ "/pp candidate events; the open histogram is 
data and the dashed line is from a ip' — > "/pp phase-space 
MC events(with arbitrary normalization), (b) Comparisons 
between data and PWA fit projection for pp mass spectrum, 
the representations of the error bars and histograms are same 
as those in Fig. [2] 



also describe the data well. For this case, the mass and 
width of the X{pp) shift by 19 MeV/c 2 and 4 MeV/c 2 , 
respectively, while the relative change in the product BR 
is 25%. These errors are considered as second (model) 
systematic errors due to the uncertainty of the model 
dependence. 

The ip' — > "/pp decay channel is also studied using event 
selection criteria similar to those used in the J/tp — > "/pp 
study. The pp mass spectrum of the surviving events 
is shown in Fig. (3][a). Besides the well known rj c and 
Xcj peaks, there is also a pp mass threshold excess rel- 
ative to PS. However, here the line shape of the mass 
spectrum in the threshold region appears to be less pro- 
nounced than that in J/ip decays. Potential background 
processes were extensively studied with an inclusive MC 
sample of 1 x 10 s ip' events and a data sample of the 
selected ip' — » ir°pp events, and these indicate that the 
pp mass threshold structure is not from any background 
source. An exclusive MC sample, generated according to 
preliminary PWA results of ip' — > ir°pp decays with BE- 
SIII data [22|], is applied to the background estimation, 
and the background level from this source in the selected 
data sample with M pp < 2.2 GeV/c 2 is determined to be 
3.4%. 

A PWA on ip' — >• "/pp which is similar to that ap- 
plied for J/ip — ¥ "/pp decays was performed to check the 
contribution of X(pp) in ip' decays and to measure the 
production ratio between J/ip and ip' radiative decays, 
R = B{iP' -> "/X(pp))/B(J/i{; ->• -fX(pp)). Due to lim- 
ited statistics of ip' events, in the PWA, the mass and 
width of X(pp) as well as its J PC were fixed to the results 
obtained from J/ip decays. Figure E^b) shows compar- 
isons between data and MC projections for the pp mass 
spectrum. As in J/ip decays, replacing the /o(2100) with 
the / 2 (2150) and the / 2 (1910) with the / 2 (1950) yields 
no significant change in fit quality. The determined prod- 
uct BR and R value are B(ip' — s- jX) x B(X ->■ pp) = 
(4.57 ± 0.36) x 1(T 6 and R = (5.08^;^)%, respectively. 

With the consideration of systematic uncertainties sim- 
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ilar to those in J/ip decays, and the uncertainty of the 
total number of ip' events, the total relative systematic er- 
ror on the BR is (tgg (syst.) ±28 (model) )% , and system- 
atic error on R values is (j^'Jg (syst.) ± 0.12 (model))%. 
Similar to all cases studied in J/ip analysis, the statis- 
tical significance of the X(pp) is larger than Q.9a in ip' 
decays. 

The PWA fits are also performed without the cor- 
rection for FSI effect. The corresponding log- likelihood 
value worsen by 25.6 than those with FSI effect in- 
cluded. The mass, width and product BR of the 
X(pp) are M = 1861 ±l(stat.) ±\ 3 (syst.) McV/c 2 , 
r = 1 ± 6(stat.) ti 8 (syst.) MeV/c 2 (a total width of 
r < 32 MeV/c 2 at the 90% C.L), B(J/ip -> 7 A(1860)) 
B(X(1860) ->• pp) = {8.6±° i (stat.)+ 2 ^ (syst.)) x 10~ 5 
and B(ip' ->■ 7 X(1860)) B(X(1860) 4 pp) = (4.15 ± 
0.39 (stat.)^'^ (syst.)) x 10~ 6 , respectively. The corre- 
sponding R value is (4.80±g;4| (stat.)^;^ (syst.))%. 

In summary, the PWA of J/ip —> jpp and ip' — > -fpp 
decays are performed. In J/ip radiative decays, the 
near-threshold enhancement X{pp) in the pp invariant 
mass is determined to be a h state. With the in- 
clusion of Julich-FSI effects, the mass, width and prod- 
uct BR for the X(pp) are measured to be: M — 
1832t,t 9 (stat.)+^ (syst.) ± 19 (model) MeV/c 2 , T = 
13 ± 39 (stat.)+i° (syst.) ± 4 (model) McV/c 2 (a to- 
tal width of T < 76 MeV/c 2 at the 90% C.L) and 
B(J/iP -> 1 X)B{X -> pp) = (9.0t?;? (stat.):^ (syst.)± 
2.3 (model)) x 10~ 5 , respectively. The produce BR 
for X(pp) in ip' decay is first measured to be B(ip' — > 



jX) x B(X -> pp) = (4.57 ± 0.36 (stat.)^;^ (syst.) ± 
1.28 (model)) x 10~ 6 and the production ratio of the 
X(pp) between J/ip and ip' radiative decays is R = 
(5.08lg;S (stat.)+°;^ (syst.) ± 0.12 (modcl))%. 

The mass of the X(pp) measured in the PWA fit with 
FSI effect included is consistent with the A(1835), but 
the width is significantly narrower. This indicates that 
either the X(pp) and the X(1835) come from differ- 
ent sources, or that interference effects in the J/ip —> 
jir + n~ri' process should not be ignored in the determi- 
nation of the X(1835) mass and width, or that there may 
be more than one resonance in the mass peak around 1.83 
GeV/c 2 in J/ip — > jir + ir~i]' decays. When more J/ip 
data are collected at BESIII, more sophisticated analyses, 
including a PWA, will be performed for the J/ip — > 77T7T?/ 
decay channel. A measurement of the relative production 
ratios for the A(1835) in J/ip and ip' radiative decays can 
further clarify on basis of their production ratios whether 
or not X(pp) and A(1835) are the same states. 
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